Primitive purine biosynthesis connects ancient geochemistry to modern metabolism.

An unresolved question in the origin and evolution of life is whether a continuous path from geochemical precursors to the majority of molecules in the biosphere can be reconstructed from modern-day biochemistry. Here we identified a feasible path by simulating the evolution of biosphere-scale metabolism, using only known biochemical reactions and models of primitive coenzymes. We find that purine synthesis constitutes a bottleneck for metabolic expansion, which can be alleviated by non-autocatalytic phosphoryl coupling agents. Early phases of the expansion are enriched with enzymes that are metal dependent and structurally symmetric, supporting models of early biochemical evolution. This expansion trajectory suggests distinct hypotheses regarding the tempo, mode and timing of metabolic pathway evolution, including a late appearance of methane metabolisms and oxygenic photosynthesis consistent with the geochemical record. The concordance between biological and geological analyses suggests that this trajectory provides a plausible evolutionary history for the vast majority of core biochemistry.

Structure and contingency determine mutational hotspots for flower color evolution.

Evolutionary genetic studies have uncovered abundant evidence for genomic hotspots of phenotypic evolution, as well as biased patterns of mutations at those loci. However, the theoretical basis for this concentration of particular types of mutations at particular loci remains largely unexplored. In addition, historical contingency is known to play a major role in evolutionary trajectories, but has not been reconciled with the existence of such hotspots. For example, do the appearance of hotspots and the fixation of different types of mutations at those loci depend on the starting state and/or on the nature and direction of selection? Here, we use a computational approach to examine these questions, focusing the anthocyanin pigmentation pathway, which has been extensively studied in the context of flower color transitions. We investigate two transitions that are common in nature, the transition from blue to purple pigmentation and from purple to red pigmentation. Both sets of simulated transitions occur with a small number of mutations at just four loci and show strikingly similar peaked shapes of evolutionary trajectories, with the mutations of the largest effect occurring early but not first. Nevertheless, the types of mutations (biochemical vs. regulatory) as well as their direction and magnitude are contingent on the particular transition. These simulated color transitions largely mirror findings from natural flower color transitions, which are known to occur via repeated changes at a few hotspot loci. Still, some types of mutations observed in our simulated color evolution are rarely observed in nature, suggesting that pleiotropic effects further limit the trajectories between color phenotypes. Overall, our results indicate that the branching structure of the pathway leads to a predictable concentration of evolutionary change at the hotspot loci, but the types of mutations at these loci and their order is contingent on the evolutionary context.

Carbon isotope evidence for the global physiology of Proterozoic cyanobacteria.

Ancestral cyanobacteria are assumed to be prominent primary producers after the Great Oxidation Event [≈2.4 to 2.0 billion years (Ga) ago], but carbon isotope fractionation by extant marine cyanobacteria (α-cyanobacteria) is inconsistent with isotopic records of carbon fixation by primary producers in the mid-Proterozoic eon (1.8 to 1.0 Ga ago). To resolve this disagreement, we quantified carbon isotope fractionation by a wild-type planktic ß-cyanobacterium (Synechococcus sp. PCC 7002), an engineered Proterozoic analog lacking a CO2-concentrating mechanism, and cyanobacterial mats. At mid-Proterozoic pH and pCO2 values, carbon isotope fractionation by the wild-type ß-cyanobacterium is fully consistent with the Proterozoic carbon isotope record, suggesting that cyanobacteria with CO2-concentrating mechanisms were apparently the major primary producers in the pelagic Proterozoic ocean, despite atmospheric CO2 levels up to 100 times modern. The selectively permeable microcompartments central to cyanobacterial CO2-concentrating mechanisms ("carboxysomes") likely emerged to shield rubisco from O2 during the Great Oxidation Event.

A productivity collapse to end Earth's Great Oxidation.

It has been hypothesized that the overall size of-or efficiency of carbon export from-the biosphere decreased at the end of the Great Oxidation Event (GOE) (ca. 2,400 to 2,050 Ma). However, the timing, tempo, and trigger for this decrease remain poorly constrained. Here we test this hypothesis by studying the isotope geochemistry of sulfate minerals from the Belcher Group, in subarctic Canada. Using insights from sulfur and barium isotope measurements, combined with radiometric ages from bracketing strata, we infer that the sulfate minerals studied here record ambient sulfate in the immediate aftermath of the GOE (ca. 2,018 Ma). These sulfate minerals captured negative triple-oxygen isotope anomalies as low as ∼ -0.8‰. Such negative values occurring shortly after the GOE require a rapid reduction in primary productivity of >80%, although even larger reductions are plausible. Given that these data imply a collapse in primary productivity rather than export efficiency, the trigger for this shift in the Earth system must reflect a change in the availability of nutrients, such as phosphorus. Cumulatively, these data highlight that Earth's GOE is a tale of feast and famine: A geologically unprecedented reduction in the size of the biosphere occurred across the end-GOE transition.

Diversity decoupled from sulfur isotope fractionation in a sulfate-reducing microbial community.

The extent of fractionation of sulfur isotopes by sulfate-reducing microbes is dictated by genomic and environmental factors. A greater understanding of species-specific fractionations may better inform interpretation of sulfur isotopes preserved in the rock record. To examine whether gene diversity influences net isotopic fractionation in situ, we assessed environmental chemistry, sulfate reduction rates, diversity of putative sulfur-metabolizing organisms by 16S rRNA and dissimilatory sulfite reductase (dsrB) gene amplicon sequencing, and net fractionation of sulfur isotopes along a sediment transect of a hypersaline Arctic spring. In situ sulfate reduction rates yielded minimum cell-specific sulfate reduction rates < 0.3 × 10-15 moles cell-1  day-1 . Neither 16S rRNA nor dsrB diversity indices correlated with relatively constant (38‰-45‰) net isotope fractionation (ε34 Ssulfide-sulfate ). Measured ε34 S values could be reproduced in a mechanistic fractionation model if 1%-2% of the microbial community (10%-60% of Deltaproteobacteria) were engaged in sulfate respiration, indicating heterogeneous respiratory activity within sulfate-reducing populations. This model indicated enzymatic kinetic diversity of Apr was more likely to correlate with sulfur fractionation than DsrB. We propose that, above a threshold Shannon diversity value of 0.8 for dsrB, the influence of the specific composition of the microbial community responsible for generating an isotope signal is overprinted by the control exerted by environmental variables on microbial physiology.

Degassing-induced fractionation of multiple sulphur isotopes unveils post-Archaean recycled oceanic crust signal in hotspot lava.

Mantle source regions feeding hotspot volcanoes likely contain recycled subducted material. Anomalous sulphur (S) isotope signatures in hotspot lavas have tied ancient surface S to this deep geological cycle, but their potential modification by shallow magmatic processes has generally been overlooked. Here we present S isotope measurements in magmatic sulphides, silicate melt inclusions and matrix glasses from the recent eruption of a hotspot volcano at El Hierro, Canary Islands, which show that degassing induces strongly negative δ34S fractionation in both silicate and sulphide melts. Our results reflect the complex interplay among redox conditions, S speciation and degassing. The isotopic fractionation is mass dependent (Δ33S = 0‰), thus lacking evidence for the recycled Archaean crust signal recently identified at other hotspot volcanoes. However, the source has an enriched signature (δ34S ~ + 3‰), which supports the presence of younger 34S-rich recycled oceanic material in the Canary Island mantle plume.

Triple oxygen isotope evidence for limited mid-Proterozoic primary productivity.

The global biosphere is commonly assumed to have been less productive before the rise of complex eukaryotic ecosystems than it is today1. However, direct evidence for this assertion is lacking. Here we present triple oxygen isotope measurements (∆17O) from sedimentary sulfates from the Sibley basin (Ontario, Canada) dated to about 1.4 billion years ago, which provide evidence for a less productive biosphere in the middle of the Proterozoic eon. We report what are, to our knowledge, the most-negative ∆17O values (down to -0.88‰) observed in sulfates, except for those from the terminal Cryogenian period2. This observation demonstrates that the mid-Proterozoic atmosphere was distinct from what persisted over approximately the past 0.5 billion years, directly reflecting a unique interplay among the atmospheric partial pressures of CO2 and O2 and the photosynthetic O2 flux at this time3. Oxygenic gross primary productivity is stoichiometrically related to the photosynthetic O2 flux to the atmosphere. Under current estimates of mid-Proterozoic atmospheric partial pressure of CO2 (2-30 times that of pre-anthropogenic levels), our modelling indicates that gross primary productivity was between about 6% and 41% of pre-anthropogenic levels if atmospheric O2 was between 0.1-1% or 1-10% of pre-anthropogenic levels, respectively. When compared to estimates of Archaean4-6 and Phanerozoic primary production7, these model solutions show that an increasingly more productive biosphere accompanied the broad secular pattern of increasing atmospheric O2 over geologic time8.

Sulfur Isotope Fractionation by Sulfate-Reducing Microbes Can Reflect Past Physiology.

Sulfur (S) isotope fractionation by sulfate-reducing microorganisms is a direct manifestation of their respiratory metabolism. This fractionation is apparent in the substrate (sulfate) and waste (sulfide) produced. The sulfate-reducing metabolism responds to variability in the local environment, with the response determined by the underlying genotype, resulting in the expression of an "isotope phenotype". Sulfur isotope phenotypes have been used as a diagnostic tool for the metabolic activity of sulfate-reducing microorganisms in the environment. Our experiments with Desulfovibrio vulgaris Hildenborough (DvH) grown in batch culture suggest that the S isotope phenotype of sulfate respiring microbes may lag environmental changes on time scales that are longer than generational. When inocula from different phases of growth are assayed under the same environmental conditions, we observed that DvH exhibited different net apparent fractionations of up to -9‰. The magnitude of fractionation was weakly correlated with physiological parameters but was strongly correlated to the age of the initial inoculum. The S isotope fractionation observed between sulfate and sulfide showed a positive correlation with respiration rate, contradicting the well-described negative dependence of fractionation on respiration rate. Quantitative modeling of S isotope fractionation shows that either a large increase (≈50×) in the abundance of sulfate adenylyl transferase (Sat) or a smaller increase in sulfate transport proteins (≈2×) is sufficient to account for the change in fractionation associated with past physiology. Temporal transcriptomic studies with DvH imply that expression of sulfate permeases doubles over the transition from early exponential to early stationary phase, lending support to the transport hypothesis proposed here. As it is apparently maintained for multiple generations (≈1-6) of subsequent growth in the assay environment, we suggest that this fractionation effect acts as a sort of isotopic "memory" of a previous physiological and environmental state. Whatever its root cause, this physiological hysteresis effect can explain variations in fractionations observed in many environments. It may also enable new insights into life at energetic limits, especially if its historical footprint extends deeper than generational.

Publisher Correction: Pelagic barite precipitation at micromolar ambient sulfate.

The original version of this Article contained an error in the barite saturation state equation in the fourth paragraph of the Introduction and incorrectly read 'Ωbarite=({134Ba2+}⋅{SO42-})/Ksp)'. The correct version removes the superscript 134 next to 'Ba2+'. This error has now been corrected in both the PDF and HTML versions of the Article.

Deconstructing the Dissimilatory Sulfate Reduction Pathway: Isotope Fractionation of a Mutant Unable of Growth on Sulfate.

The sulfur isotope record provides key insight into the history of Earth's redox conditions. A detailed understanding of the metabolisms driving this cycle, and specifically microbial sulfate reduction (MSR), is crucial for accurate paleoenvironmental reconstructions. This includes a precise knowledge of the step-specific sulfur isotope effects during MSR. In this study, we aim at resolving the cellular-level fractionation factor during dissimilatory sulfite reduction to sulfide within MSR, and use this measured isotope effect as a calibration to enhance our understanding of the biochemistry of sulfite reduction. For this, we merge measured isotope effects associated with dissimilatory sulfite reduction with a quantitative model that explicitly links net fractionation, reaction reversibility, and intracellular metabolite levels. The highly targeted experimental aspect of this study was possible by virtue of the availability of a deletion mutant strain of the model sulfate reducer Desulfovibrio vulgaris (strain Hildenborough), in which the sulfite reduction step is isolated from the rest of the metabolic pathway owing to the absence of its QmoABC complex (ΔQmo). This deletion disrupts electron flux and prevents the reduction of adenosine phosphosulfate (APS) to sulfite. When grown in open-system steady-state conditions at 10% maximum growth rate in the presence of sulfite and lactate as electron donor, sulfur isotope fractionation factors averaged -15.9‰ (1 σ = 0.4), which appeared to be statistically indistinguishable from a pure enzyme study with dissimilatory sulfite reductase. We coupled these measurements with an understanding of step-specific equilibrium and kinetic isotope effects, and furthered our mechanistic understanding of the biochemistry of sulfite uptake and ensuing reduction. Our metabolically informed isotope model identifies flavodoxin as the most likely electron carrier performing the transfer of electrons to dissimilatory sulfite reductase. This is in line with previous work on metabolic strategies adopted by sulfate reducers under different energy regimes, and has implications for our understanding of the plasticity of this metabolic pathway at the center of our interpretation of modern and palaeo-environmental records.

Electron carriers in microbial sulfate reduction inferred from experimental and environmental sulfur isotope fractionations.

Dissimilatory sulfate reduction (DSR) has been a key process influencing the global carbon cycle, atmospheric composition and climate for much of Earth's history, yet the energy metabolism of sulfate-reducing microbes remains poorly understood. Many organisms, particularly sulfate reducers, live in low-energy environments and metabolize at very low rates, requiring specific physiological adaptations. We identify one such potential adaptation-the electron carriers selected for survival under energy-limited conditions. Employing a quantitative biochemical-isotopic model, we find that the large S isotope fractionations (>55‰) observed in a wide range of natural environments and culture experiments at low respiration rates are only possible when the standard-state Gibbs free energy (ΔG'°) of all steps during DSR is more positive than -10 kJ mol-1. This implies that at low respiration rates, only electron carriers with modestly negative reduction potentials are involved, such as menaquinone, rubredoxin, rubrerythrin or some flavodoxins. Furthermore, the constraints from S isotope fractionation imply that ferredoxins with a strongly negative reduction potential cannot be the direct electron donor to S intermediates at low respiration rates. Although most sulfate reducers have the genetic potential to express a variety of electron carriers, our results suggest that a key physiological adaptation of sulfate reducers to low-energy environments is to use electron carriers with modestly negative reduction potentials.The ISME Journal advance online publication, 31 October 2017; doi:10.1038/ismej.2017.185.

Pelagic barite precipitation at micromolar ambient sulfate.

Geochemical analyses of sedimentary barites (barium sulfates) in the geological record have yielded fundamental insights into the chemistry of the Archean environment and evolutionary origin of microbial metabolisms. However, the question of how barites were able to precipitate from a contemporary ocean that contained only trace amounts of sulfate remains controversial. Here we report dissolved and particulate multi-element and barium-isotopic data from Lake Superior that evidence pelagic barite precipitation at micromolar ambient sulfate. These pelagic barites likely precipitate within particle-associated microenvironments supplied with additional barium and sulfate ions derived from heterotrophic remineralization of organic matter. If active during the Archean, pelagic precipitation and subsequent sedimentation may account for the genesis of enigmatic barite deposits. Indeed, barium-isotopic analyses of barites from the Paleoarchean Dresser Formation are consistent with a pelagic mechanism of precipitation, which altogether offers a new paradigm for interpreting the temporal occurrence of barites in the geological record.

Low viral predation pressure in cold hypersaline Arctic sediments and limits on lytic replication.

Viruses are ubiquitous drivers of microbial ecology and evolution and contribute to biogeochemical cycling. Attention to these attributes has been more substantial for marine viruses than viruses of other environments. Microscopy-based investigation of the viral communities from two cold, hypersaline Arctic springs was undertaken to explore the effects of these conditions on microbe-viral ecology. Sediments and water samples were collected along transects from each spring, from anoxic spring outlets through oxygenated downstream channels. Viral abundance, virus-microbe ratios and modelled virus-microbe contact rates were lower than comparable aqueous and sedimentary environments and most similar to deep subsurface sediments. No individual cell from either spring was visibly infected. Viruses in these springs appear to play a smaller role in controlling microbial populations through lytic activity than in marine water column or surface sedimentary environments. Relief from viral predation indicates the microbial communities are primarily controlled by nutrient limitation. The similarity of these springs to deep subsurface sediments suggests a biogeographic divide in viral replication strategy in marine sediments.

Viral Induced Microbial Mortality in Arctic Hypersaline Spring Sediments.

Viruses are a primary influence on microbial mortality in the global ocean. The impacts of viruses on their microbial hosts in low-energy environments are poorly explored and are the focus of this study. To investigate the role of viruses in mediating mortality in low-energy environments where contacts between viruses and microbes are infrequent, we conducted a set of in situ time series incubations in the outlet and channel sediments of two cold, hypersaline springs of the Canadian High Arctic. We found microbial and viral populations in dynamic equilibrium, indicating approximately equal birth and death rates for each population. In situ rates of microbial growth were low (0.5-50 × 103 cells cm-3 h-1) as were rates of viral decay (0.09-170 × 104 virions cm-3 h-1). A large fraction of the springs' viral communities (49-100%) were refractory to decay over the timescales of our experiments. Microcosms amended with lactate or acetate exhibited increased microbial growth rates (up to three-fold) indicating organic carbon as one limiting resource for the microbial communities in these environments. A substantial fraction (15-71%) of the microbial populations contained inducible proviruses that were released- occasionally in multiple pulses- over the eight monitored days following chemical induction. Our findings indicate that viruses in low-energy systems maintain low rates of production and activity, have a small but notable impact on microbial mortality (8-29% attenuation of growth) and that successful viral replication may primarily proceed by non-lethal strategies. In cold, low biomass marine systems of similar character (e.g., subsurface sediments), viruses may be a relatively minor driver of community mortality compared to less energy-limited environments such as the marine water column or surface sediments.

Sulfur isotope fractionation during the evolutionary adaptation of a sulfate-reducing bacterium.

Dissimilatory sulfate reduction is a microbial catabolic pathway that preferentially processes less massive sulfur isotopes relative to their heavier counterparts. This sulfur isotope fractionation is recorded in ancient sedimentary rocks and generally is considered to reflect a phenotypic response to environmental variations rather than to evolutionary adaptation. Modern sulfate-reducing microorganisms isolated from similar environments can exhibit a wide range of sulfur isotope fractionations, suggesting that adaptive processes influence the sulfur isotope phenotype. To date, the relationship between evolutionary adaptation and isotopic phenotypes has not been explored. We addressed this by studying the covariation of fitness, sulfur isotope fractionation, and growth characteristics in Desulfovibrio vulgaris Hildenborough in a microbial evolution experiment. After 560 generations, the mean fitness of the evolved lineages relative to the starting isogenic population had increased by ∼ 17%. After 927 generations, the mean fitness relative to the initial ancestral population had increased by ∼ 20%. Growth rate in exponential phase increased during the course of the experiment, suggesting that this was a primary influence behind the fitness increases. Consistent changes were observed within different selection intervals between fractionation and fitness. Fitness changes were associated with changes in exponential growth rate but changes in fractionation were not. Instead, they appeared to be a response to changes in the parameters that govern growth rate: yield and cell-specific sulfate respiration rate. We hypothesize that cell-specific sulfate respiration rate, in particular, provides a bridge that allows physiological controls on fractionation to cross over to the adaptive realm.

Atmospheric record in the Hadean Eon from multiple sulfur isotope measurements in Nuvvuagittuq Greenstone Belt (Nunavik, Quebec).

Mass-independent fractionation of sulfur isotopes (S-MIF) results from photochemical reactions involving short-wavelength UV light. The presence of these anomalies in Archean sediments [(4-2.5 billion years ago, (Ga)] implies that the early atmosphere was free of the appropriate UV absorbers, of which ozone is the most important in the modern atmosphere. Consequently, S-MIF is considered some of the strongest evidence for the lack of free atmospheric oxygen before 2.4 Ga. Although temporal variations in the S-MIF record are thought to depend on changes in the abundances of gas and aerosol species, our limited understanding of photochemical mechanisms complicates interpretation of the S-MIF record in terms of atmospheric composition. Multiple sulfur isotope compositions (δ(33)S, δ(34)S, and δ(36)S) of the >3.8 billion-year-old Nuvvuagittuq Greenstone Belt (Ungava peninsula) have been investigated to track the early origins of S-MIF. Anomalous S-isotope compositions (Δ(33)S up to +2.2‰) confirm a sedimentary origin of sulfide-bearing banded iron and silica-rich formations. Sharp isotopic transitions across sedimentary/igneous lithological boundaries indicate that primary surficial S-isotope compositions have been preserved despite a complicated metamorphic history. Furthermore, Nuvvuagittuq metasediments recorded coupled variations in (33)S/(32)S, (34)S/(32)S, and (36)S/(32)S that are statistically indistinguishable from those identified several times later in the Archean. The recurrence of the same S-isotope pattern at both ends of the Archean Eon is unexpected, given the complex atmospheric, geological, and biological pathways involved in producing and preserving this fractionation. It implies that, within 0.8 billion years of Earth's formation, a common mechanism for S-MIF production was established in the atmosphere.

Intracellular metabolite levels shape sulfur isotope fractionation during microbial sulfate respiration.

We present a quantitative model for sulfur isotope fractionation accompanying bacterial and archaeal dissimilatory sulfate respiration. By incorporating independently available biochemical data, the model can reproduce a large number of recent experimental fractionation measurements with only three free parameters: (i) the sulfur isotope selectivity of sulfate uptake into the cytoplasm, (ii) the ratio of reduced to oxidized electron carriers supporting the respiration pathway, and (iii) the ratio of in vitro to in vivo levels of respiratory enzyme activity. Fractionation is influenced by all steps in the dissimilatory pathway, which means that environmental sulfate and sulfide levels control sulfur isotope fractionation through the proximate influence of intracellular metabolites. Although sulfur isotope fractionation is a phenotypic trait that appears to be strain specific, we show that it converges on near-thermodynamic behavior, even at micromolar sulfate levels, as long as intracellular sulfate reduction rates are low enough (<<1 fmol H2S⋅cell(-1)⋅d(-1)).

Multiple S-isotopic evidence for episodic shoaling of anoxic water during Late Permian mass extinction.

Global fossil data show that profound biodiversity loss preceded the final catastrophe that killed nearly 90% marine species on a global scale at the end of the Permian. Many hypotheses have been proposed to explain this extinction and yet still remain greatly debated. Here, we report analyses of all four sulphur isotopes ((32)S, (33)S, (34)S and (36)S) for pyrites in sedimentary rocks from the Meishan section in South China. We observe a sulphur isotope signal (negative δ(34)S with negative Δ(33)S) that may have resulted from limitation of sulphate supply, which may be linked to a near shutdown of bioturbation during shoaling of anoxic water. These results indicate that episodic shoaling of anoxic water may have contributed to the profound biodiversity crisis before the final catastrophe. Our data suggest a prolonged deterioration of oceanic environments during the Late Permian mass extinction.

Active microbial sulfur disproportionation in the Mesoproterozoic.

The environmental expression of sulfur compound disproportionation has been placed between 640 and 1050 million years ago (Ma) and linked to increases in atmospheric oxygen. These arguments have their basis in temporal changes in the magnitude of 34S/32S fractionations between sulfate and sulfide. Here, we present a Proterozoic seawater sulfate isotope record that includes the less abundant sulfur isotope 33S. These measurements imply that sulfur compound disproportionation was an active part of the sulfur cycle by 1300 Ma and that progressive Earth surface oxygenation may have characterized the Mesoproterozoic.